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Abstract 

Monte Carlo simulations on the SK model have been done to investigate aging pro- 
cesses after a rapid quench from T = oo to the spin-glass phase. The time range of 
simulations is taken care of so that the system does not surmount free-energy barriers 
between the pure states which are considered to diverge in the thermodynamic limit. 
We simulate time evolutions of extensive energy of the system, Parisi's overlap distribu- 
tion function, auto-correlation and clones-correlation functions, distribution functions 
of the two correlations, and magnetization induced by the field applied after a certain 
waiting time. The data simulated exhibit a rich variety of aging phenomena. Most of 
them can be interpreted in a unified way, though qualitatively, by the scenario of growth 
of quasi- equilibrium domains which we will introduce in this work. More surprisingly, 
they suggest strongly that a basin of attraction of one dominant pure state spans almost 
an entire phase space of the system with a common time-reversal symmetry. 
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1 Introduction 



Since the first observation by Lundgren et al[[j]] aging phenomena in spin glasses have been 
extensively studied^. They are expected to reveal nature of the low-temperature spin-glass 
phase which has not yet been settled in spite of nearly two decades of dispute. One of 
the key concepts on aging phenomena recently introduced by Bouchaud|| and incorporated 
in the analytic theories [|Jf5| is weak- erg odicity breaking. It was originally introduced in 
the argument on the trap model ||. In a system exhibiting aging phenomena there exist 
many metastable states, or 'traps'. Depths of the traps are distributed continuously in such 
a way that the average time for the system to escape from them by thermal activation 
process is infinite; or, it takes an infinite time for the system to equilibriate. This picture 
of weak-ergodicity breaking is represented in terms of the spin auto-correlation function 
C(M') = N- 1 Ei Si(t)Si(t') as followsf|[|: 

dC(t,f) , dC(t,t') , . 

lim lim C(t,t') = q c (> 0), (lb) 

t— >oo t'— »oo 

lim C(t, f) = for a fixed t' (< oo), (lc) 

t — >oo 

where q c is a certain constant. 

Various simulations which reproduce some aspects of the experimentally observed aging 
phenomena have been already reported|J[[[]]||||. But they are not yet satisfactory to be 
able to understand the phenomena in a unified way. Since the recent phenomenological and 
theoretical arguments mentioned above, as well as recent analyses on experiments [|X0|] [|TT 
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more or less based on the mean-field picture of spin glasses, we believe it is worth examining 
aging phenomena in the mean-field spin-glass model more in detail by simulations [O]. In 
the present work therefore we have performed Monte Carlo (MC) simulation on the SK 
model [JIB]], whose equilibrium properties are well understood []TJJ. Its purpose is to obtain 



various information on its aging phenomena, including those which are related to the above- 
mentioned key concept. 

Aging phenomena we are concerned with are non-equilibrium processes toward equilib- 
rium observed in macroscopic systems in long but finite intervals of time. Since, on the 
other hand, the present simulations are performed only on finite systems with N spins, we 
have to take care of their time range. Suppose the system has a certain relaxation time, 
t erg (N,T), which may depend also on temperature T, and whose systematic dependence on 
N is known. If ]im.^^. 00 t e]: g(N, T) = oo, simulations on aging phenomena have to be done in 
the time range of t < t eTg (N, T) in order to be free from the finite-size effect. 

According to Parisi's picture[0 on the spin-glass phase of the SK model in equilibrium, 
there are many pure states which are separated from each other by insurmountable free- 
energy barriers. More explicitly, it was argued that in systems with large N there exist two 
types of the free-energy barriers: one between the pure states with a common time-reversal 
symmetry, and the other between those with opposite time-reversal symmetries ||1 5|| [ 16 ] . The 
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characteristic times, t^ rK (N,T) and t^ rg (N } T), needed to surmount respectively the former 



and latter barriers are given by (see Fig. 2 below) 

In t% g (N, T) oc iV 1/4 , and In £ g (JV, T) oc N 1/2 . (2) 

Therefore we have to set at least t max <t^ rg (N,T), where t max is the maximum observation 
time of simulations. In the present work we call such a limited time interval after a rapid 
quench aging range, and a portion of phase space that the system explores within that 
interval a 'pure state'. 

In the present work we simulate time evolution of various quantities after the quench; 
extensive energy of the system, Parisi's overlap distribution function, auto-correlation and 
clones-correlation functions|TJj[l7]][[l^], and distribution functions of the two correlations. 



The magnetization induced by the field applied after a waiting time t w is also investigated. 
It has turned out that the data simulated exhibit a rich variety of aging phenomena. They 
indicate complexity of the rugged free-energy structure within the 'pure state'. It is expected 
to have many free-energy barriers of finite heights which are surmountable within the aging 
range introduced above. Most of our simulated results can be interpreted in a unified way, 
though qualitatively, by the scenario of growth of quasi- equilibrium domains which we will 
introduce and explain in this work. Surprisingly, the results suggest that the basin of at- 
traction of one 'pure state' spans an almost entire phase space with a common time-reversal 
symmetry. We call this feature the dynamic dominant pure-states picture. On equilibrium 
properties of the SK model, on the other hand, there is the argument ]T9| that among many 



pure states only a few (even one) of them dominate the Gibbs-Boltzmann measure. We 
call this aspect the static dominant pure-states picture. Our simulated data imply that the 
dynamic dominant pure-states picture describes dynamical (or aging) processes in the phase 
space whose organization of states gives rise to the static dominant pure-states picture for 
equilibrium. 

In the next section the model and method of our simulations are briefly explained. In 
section 3 we present the results of simulations which are interpreted by our scenario of growth 
of quasi-equilibrium domains in section 4. The final section is devoted to the concluding 
remarks. 

2 Model and Method 

By means of the standard heat-bath method of the MC simulation we have studied the ±J 
SK model with mean zero and variance (N — 1) _1 / 2 [[ZIJ. The system sizes examined are 
iV = 32 ~ 2048, particularly N = 128, 512 in detail. Mainly analyzed in the present work 
are aging processes after rapid quenches from T = oo to T below the spin-glass transition 
temperature T c (= 1 in the limit N — > oo). Simulation of this process is to simply perform 
a MC run at T starting from a random spin configuration. Physical quantities are obtained 
by taking average over such M independent MC runs for each realizations of {Jij} (sample) 
and over the N s samples. Those at time t (in unit of IMC step per spin) in each MC run 
are evaluated by means of the following short time average, 

A {m \t) = -±— £ a^(t'), (3) 

^ "+* 1 t'=t-At 
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where it') is the value of quantity A at step t 1 of m-th MC run, and At -C t (eg. At = 
for t < 16 and 511 for t > 2 13 ). In the present work we put M = 10 and iV s = 40 ~ 2000 
depending on N. 

Most of our simulations are performed in the aging range mentioned in the previous 
section. Here we further make a comment on this time range in relation to the limiting 
procedures in eqs.(lb, c). For these equations the limit N —>■ oo is taken before those of 
t' and/or t. Such limiting procedures cannot be carried out in our simulational study on 
finite systems. In the present work limiting values of the types of eqs.(lb) and (lc) are 
extrapolated from the following quantities in systems with finite N: those simulated in the 
time range t «C f for the type of eq.(lb) and thoes in t ^> t! for that of eq.(lc), both t and 
t' being in the aging range. 

Another comment is on a cut-off of the aging range in the shorter side. The time range 
of simulations should not be too short in order to distinguish aging phenomena peculiar to 
spin glasses from certain initial transient processes expected to exist in any systems. We 
denote a time-scale of this short cut-off as £ q . e .. Its explicit value, as well as that of t max , will 
be discussed below referring to our simulated data. 



3 Results 

3.1 Energy and Parisi's overlap distribution 

Let us first investigate the time evolution of the energy E(t a ) after the quench at t a = 0. 
As shown in Fig. la, the energy per spin drops rapidly and then gradually saturates to the 
iV-dependent equilibrium value E eq /N. Interesting information is obtained when we look at 
how the extensive energy E(t a ) approaches to E eq . In Fig. lb we show A_E(t a ) = E(t a ) — E eq 
where E eq is approximated by E(t max ). The time t® Tg (N,T), defined by AE(t^ rg ) = T, is 
shown in Fig. 2. It has the same iV-dependence as that of t^ Tg (N,T) of eq.(2) obtained 
previously ||15||, though the proportional constants are different. In the present work we 
adopt t® g (N,T) as the upper boundary of the aging range. 

A typical barrier energy E&(t) for the system to be able to surmount within time t is 
given by £^b(^)— Tint, where the attempt time of thermal activation process is put 1 MC 
step. The value AE(t a ) + -Ee^a) for iV =512 is also shown in Fig. lb. The data shown in 
the figure tell that relaxational processes in the aging range are associated with decrease of 
AE(Q (even of AE(Q + E B (Q). 

In Fig. 3 the time evolution of Parisi's overlap distribution functions P(q; i a ) is demon- 
strated for systems with N = 128 at T = 0.4. Here q is one of MiM— 1)/2 overlaps of magne- 
tization configurations q(t) = iV -1 J2i m i k \t)mf\t) where m\ k \t) is evaluated by eq.(3) with 
a = Si. In the figure the thick broken line represents Prsb(q'), i-e., P{q) in equilibrium and 
in the thermodynamic limit which is evaluated by solving numerically Parisi's equation [pi]]. 



The position of its delta peaks are specified ±g E A; the self-overlap parameter |14]]. As already 
reported by Bhatt and Young |THj , P(q; t a ) evolves from a single gaussian centered at q = 
at t a = to a structure with double peaks (but of significant widths) similar to -Prsb(<?) 
at large enough t a . In the present analysis we determine t^ Tg (N,T) as the time at which 
P(q ~ 0;t a ) converges to -Prsb(<? — 0) and show it also in Fig. 2. We see in Fig. 3 that 
at t & = 512 ~ ^cr g (^ = 128, T = 0.4), P(g;t a ) exhibits already the double peaks around 
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q = igEAbut its peak at q ~ still remains. 
3.2 Auto-correlation function 

In Fig. 4a we show typical data of the auto-correlation function C(t + t w , t w ). It is obtained 
at T = 0.4 of systems with N = 512 in the aging range, i.e., t a = t + t w < t^ Tg (N = 
512, T = 0.4) ~ 3 x 10 4 . One can see that C(t + t w ,t w ) satisfies eq.(la), i.e., one of the 
conditions of the weak-ergodicity breaking. Also the data exhibit crossover behavior from 
relatively slow decay in t < t w to relatively rapid decay in t > t w . The behavior is common 



to C(t + t w ,t w ) simulated in various systems [|6| fl?]] || p2| , and is interpreted as a crossover 
from quasi-equilibrium to out-of-equilibrium behaviors. 

An interesting feature becomes clear if we plot C(t + t w ,t w ) against hi(t/t w ), as shown 
in Fig. 4b. As first pointed out by Baldassarri[[L2[], C(t + t w ,t w ) with different t w cross 
almost at a point t ~t crs . This behavior is certainly different from those of C(t + t w ,t w ) 



observed in the previous works |Hl|p|p2]|, whose C(t + t w , t w ) with different t w 's in the out-of- 
equilibrium range are well scaled to a single curve when plotted against t/t w . We also note 
that lnt crs ~lnt w and that the value of C(t CTS + t w ,t w ) is rather closer to 1 — T =q than ijea 
(= 0.75), where q= Jq q(x)dx is Parisi's order parameter in equilibrium [T4] . 

That the above crossing feature of C{t + t w , t w ) is a genuine property of aging processes 
in the present system is ascertained by the inspection of C(t a ,t a /2) which can be simulated 
by single aging process. As shown in Fig. 5 C(t a , t a /2) for various iV and T are nearly 
independent of t a in the range t\ < t a < t 2 (N,T). Here t 2 (^, T) is a time around which 
C(^a,^a/2) starts to decrease from its nearly constant value, and we have checked that it 
coincides with tf rg (N, T) within logarithmic accuracy. On the other hand, the shorter time 
scale t\ is regarded as t q , e , mentioned in section 2. The near constancy of C(t a ,t a /2) implies 
that C(t + t w ,t w ) with different t w 's have a nearly common value at ln(t/t w ) ~ 0. 

3.3 Clones-correlation function 



We have also investigated the clones-correlation function Q(t + t w ,t w ) |i~2[ |i~7| pq| . It is 
the correlation of two configurations starting from an identical one which has evolved up 
to t a =t w , but are evolving by means of two independent MC runs (with different sets of 
random numbers) at t a >t w . A typical result is shown in Fig. 6 for T = 0.4, iV = 512. As 
found also by Baldassarri [|l"2] |, C(t + t w ,t w ) and Q(t + t w ,t w ) cross at hit ~lnt w . Their values 
at the crossing are again closer to q independently of t w . A more remarkable feature seen 
in the figure is that Q(t + t w ,t w ) with fixed t w tend to saturate to constant values which 
depend on t w as t — > £c rg (iV, T). It is also noted (not shown) that Q(t + t w , t w ) plotted against 
ln(t/t w ) also cross nearly at a point similarly to C(t + t w ,t w ) in Fig. 4b. 

3.4 Fluctuation-dissipation theorem 

In order to get further insights into aging processes we have simulated the magnetization 
m(t; t w ) induced by the field h which is switched on at t a = t w . The result for T = 0.4, iV = 
512 and with h = 0.1 is drawn against ln(t/t w ) in Fig. 7. Also shown by broken lines in the 
figure are m c (t; t w ) defined by 

m c {t;t w ) = h{l-C{t + t w ,t w )}/T. (4) 
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The fluctuation-dissipation theorem (FDT) tells that, in equilibrium and for an infinitesi- 
mally small h, m(t; t w ) and m c (t; t w ) are independent of t w and m{t) = mc(t) holds. We see 
in Fig. 7 that, as has been also observed in other model systems0 P[P2], m(t; t w )=m c (t; t w ) 
holds in the time range t <Ct w with t w larger than a certain value. The latter value, which 
roughly coincides with t\ introduced above, is also used to estimate t q . e . which specifies the 
initial transient regime of the aging process. 



4 Discussions 

4.1 Scenario of growth of quasi-equilibrium domains 

In order to interpret the above results of our simulations, let us first remind of the following 
properties of the SK model in equilibrium. It is known that the TAP equation of the 



SK model has a huge number of solutions, or local free-energy minima, in the spin-glass 
phase [25). In equilibrium, however, only a very limited number of the local minima have 



significant probability weights P a (oc exp(— F a /T)) where F a is the (extensive) TAP free- 
energy of the a-th solution |19| (the static dominant pure-states picture mentioned in §1). 
Associated with such lowest free-energy states are limited portions of phase space which are 
separated by free-energy barriers described by eq.(2). One of the possible interpretations is 
then that aging processes simulated in the present work reflect the free-energy structures of 
such local regions centered at each lowest free-energy state. We may identify these regions 
to the 'pure states' introduced in section 1. 

Based on the above assumption we argue that each 'pure state' contains a huge number of 
local minima (which we call hereafter simply states), as schematically shown in Fig. 8. The 
lower their (extensive) energies, the less is the number of the states. [Here and hereafter we 
disregard entropy effects which have not been evaluated explicitly in the present simulation.] 
Their basins of attraction for a rapid quench process are considered to have more or less 



comparable sizes ||26|| . Therefore after a certain initial transient time the system is found, 
with probability of almost unity, in one of the huge number of states with a relatively higher 
energy, from which it starts to look for lower and lower energy states by thermal activation 
processes; i.e., it ages. We consider that t q , e , introduced before corresponds to the above- 
mentioned transient time, it is noted that t q . e . are little dependent on N and is less than a 
few tens MC steps (see Figs. 5 and 7). 

Let us then suppose that at time t & = £ w i in the aging range the system reaches to a 
certain state, say Si in Fig. 8. In a time interval of t & = t w i + 1 with t <C i wl the system 
is fluctuating in a local region Ri centered at Si (the shaded area in Fig. 8), thereby it is 
expected to visit states in the region with frequencies according to their relative Boltzmann 
weights. We call such regions quasi-equilibrium domains. As time goes on another interval of 
order £ wl , i.e., at t a = t w2 — 2t wl , the system surmounts a barrier of order Tlnt wl and finds a 
lower energy state, say S 2 in Fig. 8. This corresponds to entering into the out-of-equilibrium 
range of the waiting time i wi . The state 5 2 is, in turn, associated with the quasi-equilibrium 
domain R2 of the waiting time t W 2 which is larger than Ri. The quasi-equilibrium domain 
grows with time until it exhausts the whole 'pure state', or the system reaches the lowest 
energy state Soo in Fig. 8 (which takes an infinite time if N = 00). 

The above scenario for aging processes, which we call growth of quasi-equilibrium domains, 
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has been initially introduced from our relaxational-modes analysis on the same problem [27]. 
From the present MC study it is introduced to explain straightforwardly the following results 
in the aging range of t q e . < £ a < t^ rg (N, T): 1) the monotonic decrease of the extensive energy 
in the aging range (AE(t a ) in Fig. lb), 2) the crossover from quasi-equilibrium to out-of- 
equilibrium behaviors at around t ~i w seen in C(t + t w ,t w ) of Fig. 4, 3) the apparent FDT 
at t <Ct w demonstrated in Fig. 7, and 4) the behavior of the clones-correlation shown in 
Fig. 6 which is explained as follows. Suppose the clones are created at t a =t wl from a spin 
configuration near S± in Fig. 8. Then at t > £ w i, the clones visit more frequently near 
the same state S2 whose relative Boltzmann weight is largest in the larger domain R 2 . The 
Hamming distance between them is then considered to be smaller than that between Si and 
5*2 which determines predominately the value of the corresponding auto-correlation function. 
This means Q(t + £ w i, t w i) > C(t + t wi , t w \) at t > t w i as observed in the simulation. 

Concerned with the above items 2) and 3), we note that at t <Ct w C(t + t w , t w ) depend not 
only on t but also on t w (see Fig. 4a). This means that nature of overall fluctuation within 
each quasi-equilibrium domain differs between the domains and from true equilibrium. Still 
C{t + t w , t w ) and m(t; t w ) obey the apparent FDT. This feature is certainly different from the 
recent arguments f5f[ [|T0| that both C(t + t w , t w ) and m(t; t w ) consist of two parts, equilibrium 
and aging ones, and that the FDT holds only for the equilibrium parts which depend on t 
alone. 

We also note that the above time range t w ^> t( > t q . e .) is the one where we look for the 
limiting value of eq.(lb) as pointed out in Section 2. Naive speculation based on the results 
shown in Figs. 4b and 5 is that q c in eq.(lb) is certainly less than qea i n contrast to the 
previous claim of ? c = ?ea[0- In our simulation on a large but finite N system, we regard that 
the system is in 'equilibrium' within one 'pure state' at t w of the order of (N, T). The 
extrapolation to iV = 00 keeping these conditions at t w ~tg rg (iV, T) gives rise to q c smaller 
than q-EA- 

4.2 Distribution functions of clone- and auto- correlations 

The time evolution of distribution functions P(Q;t,t w ) with Q(t) = N^ 1 J2i m ^{t) m t\t) 
shown in Fig. 9a further supports the above interpretation of the clones-correlation. Here 
mf\t) is evaluated by eq.(3) whose a in the r.h.s. is Si of the k-th clone. The distribution 
is therefore obtained from MN S data points of Q{t). A sharp peak of P(Q;t,t w ) at t ^ t w 
located around Q ~ ^ea indicates simply that the clones are not yet separated enough and 
are fluctuating near the lowest energy state in the quasi-equilibrium domain of waiting time 
t w . Around t ~ t w height of the peak decreases and its width increases. This is interpreted 
that the clones now tend to escape statistically independently of each other from the domain 
of £ w . Interestingly at t ^> t w the peak near Q ~ qea sharpens again though its height does 
not recover very much. The peak in this time range suggests that the clones are again near 
the same lowest energy state but of a larger quasi-equilibrium domain than the original one. 
This implies that Q(t + t w ,t w ) for fixed t w converge to certain finite values. This is in fact 
the case as we have already seen in Fig. 6; Q(t + £„, t w ) almost saturate to certain constants. 
Correspondingly, the shape of P(Q;t,t w ) becomes little dependent on t at t (~tl s (N,T)). 

In order to examine the iV-dependence of the above characteristic feature in the clones- 
correlation, we show in Fig. 9b Q(t + t w , £ w ) at t (N, T) for N =128 and 512 both with 
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t w =64. The two almost coincide with each other. The small difference in the peak positions 
near Q =c/ea is attributed to the finite-size effect. Naive extrapolation of these results to 
iV — > oo (thereby t^. (N, T) — > oo) is lim t __ >00 (5(t + t w , £ w )> for a fixed t w (< oo). Thus the 
aging process within one 'pure state' of present interest belongs to type I by means of the 
classification introduced by Barrat et al|jTB|. 

In Fig. 9b we also show, by the thin solid curve, Q(t+t w , t w ) at t ~tg rg (iV, T) for N = 128. 
One sees that a small peak starts to develop not around Q = but around Q = — Qea- This 
result is interpreted as follow. When one of the clones surmounts the energy barrier between 
the phase spaces with different time-reversal symmetry, it goes down and situates near the 
counterpart of within a time interval of t^. (N, T)(<^tg rg (iV, T)). 

In contrast to Q(t + t w ,t w ) above discussed, C(t + £„,£„) in Fig. 6 doesn't exhibit 
tendency of the saturation to certain finite values within the aging range of the present 
concern. Although it neither completely vanishes in the same time range, we expect that 
the weak-ergodicity breaking condition of eq.(lc) holds in the limit N oo within its 
aging range. This expection is derived from inspection of the distribution functions of auto- 
correlation, P(C;t,t w ) with C(t) = N^ 1 Ei Si\t + t w )Si(t w ), shown in Fig. 10. At t > t w , 
in contrast to P(Q] t, t w ), P{C\ t, £ w ) tends to become a gaussian form whose center tends to 
approach C — 0. We therefore consider that the weak-ergodicity breaking of eq.(lc) holds 
since any configuration at t a =t w (< oo) in a 'pure state' is orthogonal to (separated far 
away from) the configuration at t a = oo, i.e., the lowest energy state of the 'pure state'. 

According to Barrat et al|18[, the above aging process of type I judging from Q(t + t w , t w ) 
and with the weak-ergodicity breaking judging from C(t + t w , t w ), is expected to appear, for 
example, in a system having a 'gutter' in phase space. But our scenario of growth of quasi- 
equilibrium domains is also compatible with the two characteristics. We only need a quite 
natural assumption that a basin of attraction of a 'pure state' is infinitely large if N = oo. In 
fact our scenario is more appropriate judging from the time evolution of P(Q;t,t w ) already 
described above. We may say that the energy structure in one 'pure state' is, instead of a 
'gutter', a 'funnel' with an infinitely wide input-mouth. 

4.3 Dominant pure-states pictures 

Lastly let us compare the three distribution functions P(Q',t ~ t^ rg ,t w ), P(q;t a ~ tf rg ) and 
Prsb{q) drawn in Fig. 9c. Here P(q;t a ~ L) and -Prsb(<?) (also those in Figs. 9a,b) 
are normalized as f_ x P(q)dq = 2, differently from those in Fig. 3. The reason of this 
normalization is that time evolution of the clones at t ^t^ Tg (N,T)(<^.t^ Tg (N,T)) of present 
interest is almost confined in a part of phase space with a common time-reversal symmetry. 
[A small weight of P(Q;t,t w ) smearing out to Q < is considered due to the finite-size 
effect]. It is rather surprising that near coincidence, P(Q; t ~ t^ rg , t w ) ~ P(q; t a ~ ^(A^, T)) 
— -Prsb(<?), holds for Q,q > 0.3. [To be more accurate, near equality P(q; t a ) = Prsb(?) for 
all q is obtained only at t larger than t^ rg (N, T) as seen in Fig. 3 and in the previous work[[il|]. 
By 'near coincidence' here we disregard quantitative difference in the weights associated with 
the peak at Q, q ~^ea and qualitative difference in the shapes at Q, q < O.3.] 

The near coincidence P(Q;t ~ t^ rg ,t w ) ~ Prsb(q) alone indicates that the organization 
of states in one 'pure state' centered at the lowest energy state (S^ in Fig. 8) is similar to 
that of the whole phase space with a common time-reversal symmetry which gives rise to 
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Prsb(q)- The other near coincidence P(Q',t ~ t^ rg ,t w ) ~ P(q;t& ~ tf rg ) tells more since its 
r.h.s. quantity is the overlap between two configurations which start from random initial 
configurations independently chosen. Its plausible interpretation we can think of is that a 
basin of attraction of one 'pure state' reached by our simulation in fact covers almost an 
entire phase space with a common time-reversal symmetry. This is what we have called 
the dynamic dominant pure-states picture in section 1. We regard the near coincidence of 
the three distribution functions as an evidence that the aging processes we have simulated 
are stochastic dynamics of the system in the phase space, whose organization of states in 
equilibrium is described by the static dominant pure-states picture|T^]. In other words, our 
present simulations yield, for the first time to our knowledge, an evidence for the latter 
picture from a dynamical point of view. 



5 Conclusion 

We have proposed the scenario that aging processes at a fixed T (and h) in the SK model, 
which are observed in the aging range of t q . c . < t < (N,T), are growth processes of 
quasi-equilibrium domains, or stochastic dynamics of the system looking for dominant pure 
states which exhaust the probability weights in equilibrium. It qualitatively explains most 
of our simulated results. In our arguments extensive energies of states (solutions of the 
TAP equation), which are represented by the ordinate of Fig. 8, and which govern relative 
Boltzmann weights in each quasi-equilibrium domain, play a central role. If, however, we 
want to know more details of aging phenomena, such as an explicit functional form of C(t + 
t w ,t w ), we have to analyze the organization of the states in the direction of the abscissa of 
Fig. 8. For this purpose the difference in the extensive energies may not be so important than 
the above argument (imagine that in the SK model with iV — > oo the system is somewhere 
in the middle of the 'funnel' of Fig. 8 at a finite time). The research in this direction is now 
going on. 
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Figure Captions 



Fig.l: Time evolution of energy at T = 0.4. a) Energy per spin in systems with N(N S ) =32(2000), 
64(2000), 96(2000), 128(1000), 192(600), 256(500), 512(200), 1024(40), and 2048(50) 
from top to bottom, b) Extensive energy relative to E eq for iV = 32 ~ 512. The data 
with broken line represents AE(t a ) + E^tg) for iV = 512 (see text for definition of 
E B )- 

Fig.2: The system-size (N) dependence of characteristic relaxation times related to aging 
processes in the SK model at T = 0.4 for N < 256; t^ g (JV,T) and t s eig (N,T) of eq.(2) 
and t% JN,T) determined by AE(t® ) = 0.5 (= 1.2T) from Fig. 2b. The solid line is 



ln(tg (AT, T)) = 2.56N 1 / 4 - 0.66 due to Ref.pl . The area below JN, T) is the aging 



range introduced and studied in the present work. 

Fig.3: Parisi's overlap distribution function P(q; t a ) for A" = 128, T = 0.4. It is normalized 
as /_! P(q; t)dq = 1. The data are t a =16, 64, 512, 4096 and 32768 from top to bottom 
at around q ~ 0. The thick broken line represents -Prsb(<z)- 

Fig. 4: Auto-correlation functions C(t + t w ,t w ) with t w = 2 n (n = 3 ~ 13) plotted against t (a) 
and t/t w (b) (T =0.4, N =512, N s =200). 

Fig.5: Auto-correlation functions C(t a , t a /2) for some T and N. 

Fig.6: Clones-correlation functions Q(t + t w ,t w ) with t w =8, 16, 64, 256 and 1024 (T =0.4, 
A^ =512, N s =100). Auto-correlation functions C(t + t w ,t w ) are also shown by lines 
with the open symbols. The arrow indicates t =t^ Tg (N,T). 

Fig. 7: Induced magnetization m(t; t w ) (symbols with solid lines) and mc(t; t w ) (broken lines) 
of eq.(4) with t w = 2 n (n = 3 ~ 13) plotted against t/t w (T =0.4, h =0.1, A^ =512, 
A^ s =200). 

Fig.8: Schematic representation of the rugged energy structure of one 'pure state' of the SK 
model. The dots represent local energy minima (Si, S 2 , ■ ■ ■ , Soo) and the shaded areas 
quasi-equilibrium domains (Rx, R2, • • •)■ 

Fig.9: a) Distribution functions of clones-correlation P(Q;t,t w ) corresponding to the curve 
with t w =64 in Fig. 6. The broken curve without symbols represents -Prsb(<?) (also in 
Figs. 9b,c, and 10) b) Comparison of P(Q;t,t w ) with N = 128,512 and t w =64: the 
solid curve with circles represents P(Q] t = 512 ~ t® Ig , t w ) oi N = 128, the broken 
curve with triangles P(Q;t = 2 15 ~ t^ rg ,t w ) of N = 512, and the thin solid curve 
P(Q; t = 2 15 ~ t^ rg , t w ) ofN = 128. c) Comparison of P(Q; t ~ tf rg , * w ) (broken curve 
with squares), P(q;t a ~ t^ rg ) (solid curve), and ~ Prsb(q)- f° r AT = 512 and T = 0.4. 
In this figure P(q) are normalized as j\ P(q)dq = 2. 

Fig. 10: Distribution functions of auto-correlation P(C;t, t w ) corresponding to the curve with 
t w =64 in Fie. 6. 
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